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binodal for a given interaction potential of Lennard- Jones type. We use this correspondence 
further to construct a distribution of the Bose-Einstein type for a classical gas with the 
help of the new notion of Bose condensate, making it possible to decrease fractal dimension 
while simultaneously preserving the number of particles. In so doing, we use new global 

i-^ . distributions in number theory. 
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This paper is dedicated to Stephen Smale on the occasion of his eightieth birthday. 

1 Introduction 



(^ ' At the origin of the mathematical theory of distributions (generahzed functions) which result 

from the fundamental works of L. Schwartz, and also of Gel'fand and Shilov, the so-called 
"Schrodinger problem", or the "Schrodinger conjecture", arose in mathematics. In the famous 
KJi \ memoir in which Schrodinger introduced his equation, he gave definitions of eigenfunctions for 

^ ■ the discrete and the continuous spectrum. He defined the eigenfunctions of the discrete spectrum 



correctly, whereas the definition of functions corresponding to the continuous spectrum contained 
an inessential error. Namely, Schrodinger assumed that the functions corresponding to the points 
of continuous spectrum are bounded at infinity. However, in fact, these functions can have a 
growth which cannot exceed that of some power of r, where r stands for the radius, as r — )• cxd. 

I published the related counterexample in 1968 |15) . The mathematical proofs of the above 
conjecture (such proofs were presented by Maurin, Kostyuchenko, Gel'fand, and Shilov) con- 
tained errors, which the authors had found themselves when studying my counterexample. A 
counterexample in the case of an absolutely continuous spectrum was given in 1993 by the author 
of the present paper and S. Molchanov and was reported in a plenary talk at the Congress on 
Mathematical Physics. The so-called Steklov problem for polynomials was solved simultaneously 
(in the negative; see [35]). 

Certainly, I had no doubt that Schrodinger's definition of eigenfunctions is inexact; however, 
I presented the counterexample only after the problem in question passed to the area of math- 
ematics. The point is that there are many well-known paradoxes in physics. On the contrary, 
there must be no paradoxes in mathematics (in mathematical physics). The well-known Gibbs 
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paradox was regarded as a physical paradox, and fifteen Nobel prize laureates tried to solve 
it. However, two famous mathematicians, von Neumann and Poincare, also tried to solve this 
paradox. After the exhausting paper by V. V. Kozlov [9j, the Gibbs paradox finally moved to 
the mathematical area. Kozlov proved the following fact [TO]: it follows from the Poincare model 
concerning the dynamics of collisionless gas in a rectangular parallelepiped with mirror walls 
that, if an interior wall between two parts of the vessel disappears, then the entropy increases 
stepwise. 

However, it is clear that, if we evaluate the density, i.e., the number of particles in a unit 
volume, then the entropy, both as the logarithm of the number of possible variants and as the 
Kolmogorov complexity, is preserved. 

The problem of correctly defining the notion of ideal gas became now mature in mathematics 
as well. 

It should be noted that the microcanonical Gibbs distribution holds (see Theorem 2). 

A completely different approach not connected with the microcanonical Gibbs distribution 
was proposed simultaneously by Green, Kirkwood, and others. This approach was studied in 
great detail by N. N. Bogolyubov. However, in his construction, Bogolyubov used the conjecture 
of chaos preservation stated by Mark Kac. A counterexample to this conjecture was constructed 
by the author of the present paper and O. Shvedov in |40j . 

2 Partition of integers and revision of the Bose— Einstein distri- 
bution 

The relationship between the Bose-Einstein distribution and number theory was studied, in 
particular, in [1], [36], [37], [38], [12], [l3]. Let us consider two examples. 

As an example of a simple model of Bose condensate, consider Koroviev's trick (well known 
from M. Bulgakov's novel "Master and Margarita" ) of scattering money bills in a variety theater's 
audience. According to number theory and Kolmogorov complexity, if Koroviev had one million 
bills and the number of spectators were ten thousand, then only one thousand of them would get 
bills (see [37]). The other nine thousand would not get a single bill (and presumably would die 
of hunger). This is exactly a model of Bose condensate. But if the spectators united into groups 
of ten and agreed to divide the bills between themselves, then, figuratively speaking, no one 
would die, i.e., the number of spectators (or particles) would be preserved. But the association 
of the spectators into groups would mean a constraint on their degrees of freedom, just as the 
association of particles into clusters. 

This leads to two conclusions. 

First, from the mathematical point of view, this example is equivalent to the existence of 
a two-dimensional Bose condensate, but this refutes the physicists' postulate (dating back to 
Einstein) that no two-dimensional condensate does exists. It turns out that, in the Bose-Einstein 
distribution, it is necessary to add a special term taking into account the fact that the number of 
particles is finite. Simultaneously, it provides the asymptotic distribution function (previously 
not known) for the number of bills obtained by groups of spectators. 

Second, the Bose condensate can be regarded as an association of dimers, trimers, and 
clusters, not only of "frozen" particles precipitating to the Bose condensate. And this applies 
to a classical gas, not a quantum, one. Thus, under the condition N = const, the Bose-Einstein 
distribution can be regarded as a distribution for a classical gas and the degeneracy temperature 
as the critical temperature [23]. 

Example 1. By way of an example, similar to approach of molecular dynamics, let us consider 
a set of unnumbered billiard balls of unit mass and the same color. First, consider one billiard 
ball and launch it from some (arbitrary) point with velocity not exceeding a certain sufficiently 
large value v, i.e., with energy not exceeding v'^/2. However, since the computer has certain 
accuracy, it follows that the energy of the ball will take a finite integer number s of values in 



the interval [0, i;^/2]: Aj = iEq, i = 1, . . . , s, of energies where Eq corresponds to this accuracy. 
Thus, we obtain the spectrum of energy values which can be regarded as a self-adjoint diagonal 
matrix of order s, where s 2> 1. 

By assigning such a discrete set of energies to a ball, we obtain the wave-particle correspon- 
dence in classical mechanics, because the resulting matrix is unitarily equivalent to any operator 
L with such a spectrum in a Hilbert space H. 

How many balls must be launched so that the computer is not able to determine their initial 
data? 

The spectrum corresponding to N balls can be obtained by considering the tensor product 
of N Hilbert spaces H and the corresponding spectrum of operator 

Lm = L®L®---®L (iV times). 

The eigenvalues of this operator are of the form 



£ = Y^NiXi 



i=l 

If we only consider the eigenfunctions, symmetric with respect to the permutation of the par- 
ticles, of this operator, which corresponds to the identity of the balls, then the eigenvalue 
£ = ^l^i NiXi is of multiplicity equal to the number of all possible variants of the solution of 
the problem 

s s 

£ = Y,NiXi, ^iVi = Ar. (1) 

If d is the topological dimension, then 

_ (z + d-1)! 
^' - i\id-l)\ ■ ^^) 

The binary logarithm of solutions of these Diophantine equations is called Hartley entropy. Let 
us consider a problem 

s s 

J2NiXi<£, J2Ni = N. (3) 

j=l i=l 

The entropy of problems ([T|) and (^ coincide up to \/Nln N (the accuracy up to which we solve 
these problems). 

This consideration allows to extrapolate the above theory of integers to non-integral dimen- 
sions. We shall consider relations of the form 

y- r(^ + ^) N<£ (4) 

oo 

^N, = N, (5) 



where T{d) is the F-function, and the number of solutions satisfying inequality (JU and equal- 
ity ((H) for non-integer d ("fractial dimension"). 

Since the initial set of energies is "without preferences", i.e., is in a general position, then 
all the multiplicities corresponding to ([T]) are equiprobable. The computer calculation time is 
related to the computer accuracy Eq with respect to energy. The problem is how to determine 
the number N for which, at a particular instant of time, the computer cannot recover the initial 
data in view of inaccuracy of the classical pattern or that of the quantum mechanical pattern 
(which is more accurate, but more cumbersome); the latter pattern, in turn, is not accurately 
described by an interaction of Lennard-Jones type. 



Thus, we can draw the following conclusion. The initial data in classical and quantum 
mechanical problems are discarded due to external noise. As a result, the problem is reduced 
to the distribution of {Ni} in ([T|). In this problem, we assume a priori that the initial data 
are discarded and, therefore, so is the numbering of classical particles. Although the textbook 
[TT falsely interpreted the difference between the quantum mechanical pattern and the classical 
onqj, it, nevertheless, gave a valid interpretation of the numbering of identical balls. Therefore, 
we can take symmetric eigenfunctions for Ljv- 

Therefore, it only remains to obtain the distribution of number Nj of particles using rela- 
tions ([l]). If s 3> iV, then relation ([l]) can be expressed as 

oo oo 

EoJ2^Ni = £, Y,Ni = N- (6) 

These relations coincide with those in the classical number-theoretic problem under the condition 
that £/Eq is an integer, which, of course, is of no importance in the asymptotics in s — t- oo and 
iV-^ oo. 

Thus, since the noise component has prevented us to recover the initial datao and the number 
of particles is preserved and so is the total energy £ (at least, the latter is not increased), without 
giving any preference we assume all the variants satisfying the relation 

oo J, oo 

to be equiprobable. 

Under this approach, we can also take into account collision of billiard balls. Indeed, the 
initial energy of all balls can only decrease due to friction and the passage of kinetic energy into 
thermal energy during collisions. The number of balls will remain the same and the total energy 
will not exceed the initial energy £. 

There is a similar problem in number theory. Let n be a positive integer. By a partition of 
n we mean a way to represent a natural number n as a sum of natural numbers. Let p{n) be 
the total number of partitions of n, where the order of the summands is not taken into account, 
i.e., partitions that differ only in the order of summands are assumed to be the same. The 
number Pk{n) of partitions of a positive integer n into k positive integer summands is one of the 
fundamental objects of investigation in number theory (see [36], [42], |43j ) . 

For a given partition, denote a number of summands (in the sum) equal to 1 by A^i, the 
number of summands equal to 2 by A'2, etc., and the number of summands equal to i by A'^j. 
Then X^ Aj = /c is the number of summands, and the sum Yl ^^i is obviously equal to the 
partitioned positive integer. Thus, we have 

00 00 

Y,^Ni = n, Y,Ni = k, (8) 

where Aj are natural numbers not exceeding k. 

These formulas can readily be verified for the above example. Here, all families {Aj} are 
equiprobable. 



^The authors explain the identity principle for particles as follows: "In classical mechanics, identical particles 
(such as electrons) do not lose their 'identity' despite the identity of their physical properties. ... we can 'number' 
them and then observe the motion of each of them along its trajectory; hence, at any instant of time, the particles 
can be identified ... In quantum mechanics, it is not possible, in principle, to observe each of the identical 
particles and thus distinguish them. We can say that, in quantum mechanics, identical particles completely lose 
their 'identity' [11], p. 252. (But, as a matter of fact, if the initial data for the Cauchy problem does not possess 
a symmetry property, then the situation in quantum mechanics does not differ from that in classical mechanics). 
As to the well-known discussion between Boltzmann and some mathematics [14], of course, if the particles of 
the gas are distinguishable and can be numbered, then they can also be turned back and returned to their initial 
state. 



A new global distribution for Ni <k and ^N-i = k is determined from the relations 

Vf \ ^ \=k VT '- - \=n (9) 

where 6 > and k > are constants defined from (l9|), n/k is sufficiently large, numbers n and 
k are also large, and, by using the Euler-Maclaurin summation formula, we can pass to the 
integrals (for the estimates of this passage, see [39]): 

oo 

1 ^ > , , 

drc ^ fc, (10) 



gb(x+K) ^ ^bk{x+K) -^ 





' dx^n. (11) 



g6(x+K) 2 Qbk{x+K) -^ 



It can be proved that k = gives number /cq with satisfactory accuracy. Hence, 

oo 


Consider the integral (with the same integrand as in (|12|) ) taken from e to oo and perform 
the change of variables te = .^ in the first term and hkQX = ^ in the second term. Passing to the 
limit as e — )■ 0, we derive: 

oo oo \ ebk(t 

kf) = — lini I / —^ : — / —^ r I = — lini 



b e^+0 \ J e^ — 1 J e^ — 1 j b e^+o J e? — 1 

,£& ebko / eb 

ebko 



1 r ^ic 1 1 

- lim I ^ = T hm \hi{ebkn) — hi{eb)\ = -lukn. (13) 

b £^+0 J ^ b £^+0 L V u; V ;j ^ u K / 



eb 

On the other hand, making the change 6x = ^ in (lllj) . we have: 

1 r ^d^ ^ 



62 7o e? - 1 
Relations ([H]) and ([M]) imply 



(14) 



"'ini ^V ^■o = ^ln..(l+o(l)). (15) 

Let us now find the next term of the asymptotics by setting 

ko = c^ n ' Inc^ n ' + an ' + o{n ' ), where c 



V6- 



Using the formula 

ko = c~ n ' In /cq 

and expanding In/co in 

a 



c ^ Inc ^n^/2 
we derive 

a = -21n^. 



Thus, we have obtained the Erdos formula [3]. 

We note that, in the Koroviev's trick mentioned above, the problem is put not quite similar 
to that in number theory on partitions of a positive integer n into k positive integer summands. 
Namely, the problem is put in the following way - for a given positive integer number k, find 
k < k for which the number pk{n) of partition variants is maximal. We have: (i) if fc < feo, then 
the maximal number of variants corresponds to k = k; (ii) ii k > kQ, then the maximal number 
of variants is achieved for k = k^. This is precisely the Bose-condensate (see [6j, |28] . [41j). 
Example 2. The above case corresponds to the consideration of the topological dimension 
d = 2. Consider now the one-dimensional case of a Bose condensate, which is important in 
physical problems. In the notation used in statistical physics, £ is the energy, ie are energy 
levels, k corresponds to a number N of particles and n corresponds to £/e. 

Define constants b and k from the following relations: 

oo 

^ ( eM^+!)-l - ^m^^^ dVe = J = n, (16) 











1 N \ r- 

dy^ = N. (17) 



oK^+K.) _ 1 f.bN{^+K) _ I 



Just as in Example 1, we take n, N and n/N to be sufficiently large, Ncr = ko. 
First, since 

oo oo 

_ 1 f Vldi 1 f Vldi 



2j e*«-l 263/2y e«-l' 



we have: 

For Ncr (corresponding to k = 0), we obtain the following relation: 

^ Ncr 

obS. _ I ^bNcri _ \ 




(18) 



(19) 



^cr=l[-r^--jj^]d^i (20) 







1 /■/ 1 Ncr 



VbJ Ve«'-1 e?'^- 





dc 



Subtracting 1/^ from both terms in (pO 



1 



Ncr = — I[ ^r—r - 7^ Me+ — / I 7^ - .... . ;v...o. 1 dC 




and using relations 




Ncr-eii + ^e) 



iV. 



dt (21) 



Ncr^ii + ^e) 



d^ (22) 



.2 , 



vHi + \ 



dr], 



1 



(23) 



7?2(l + f) r]' 2(l + f) 
we come to the quadratic equation 

iy^rf-Wy^r + W = (24) 



with respect to \/Ncr, where 



-1/3 



Solving this equation, we obtain: 




This imphes that, for large n, 
where 



(26) 



Ncr « c^n^/s (27) 



r ?5,-3^ < 



1 r°° V^di \ ^ 

2 Jo 7^^ J 

We stress that A'^cr determines the transition to the Bose condensate and plays the same role 
as ko (see (J12p ) in number theory. This problem in number theory is connected with the Waring 
problem. 



* * * 



In statistical calculations of the number of inhabitants in a town, the permutation between 
a child and an old man does not change the total number of inhabitants. Hence, from the point 
of view of the statistics of the given calculation, they are indistinguishable. From the point of 
view of the experimenter who observes the molecules of a homogeneous gas using an atomic 
microscope, they are indistinguishable. He counts the number of molecules (monomers) and, 
for example, of dimers in a given volume. Dimers constitute 7% in the total volume of the 
gas (according to J.M.Calo). This means that the experimenter does not distinguish individual 
monomers just as dimers and counts their separate numbers. His answer does not depend on 
the method of numbering the molecules. 

These obvious considerations are given for the benefit of those physicists who relate the fact 
that quantum particles are indistinguishable with the impossibility of knowing the world. I do 
not intend to argue with this philosophical fact, but wish to dwell only on mathematics and 
statistics and distributions related to the number of objects (compare [2], [22]). 

The existence of the A-point in Helium-4 is not a quantum effect. According to the ex- 
perimental data, the compressibility factor Z decreases to Z = 0.00763, and hence the fractal 
dimension decreases to almost 2 (see below). The derivative of energy with respect to T contains 
the derivative with respect to the chemical potential as a summand. 

This term yields a logarithmic singularity as iV — t- oo (just as in |28j). Moreover, if the 
pressure is equal to approximately 0.05 atm., then there is no singularity in the density of the 
gas. 

An example of a nonholonomic constraint (cluster) is a sphere containing another sphere 
gliding on two parallel-oriented skates fastened to its diameter. Such a collection of gliding 



spheres inscribed in each other is an example of a nonholonomic cluster. Here the larger sphere 
must coincide with all the spherical volume of a drop (see below), and the number of spheres 
must correspond to the number of constraints. As temperature decreases below T = 2.172, there 
will be an increase in nonholonomic constraint and decrease in the number of degrees of freedom 
according to the first parameter. 

Remark 1 In author's papers [IT], [E], [19] it is shown that superfluidity is not a pure 
quantum effect, but a semiclassical one. 

By the Feynman scheme (as is well known, he treated a positron as an electron moving 
back with respect to time), a particle which reflected when moving back with respect to time 
collided at the boundary with a particle which moved directly with respect to time. These two 
particles annihilated and gave a part of energy to the boundary. One can readily imagine what 
does it mean after processing some fllm inversely (back with respect to time). The tangential 
component, under this understanding of reflection, corresponds to the flow of the "pair" along the 
boundary. Therefore, to take the boundary of the vessel into account (especially if we consider 
a gas in which the convection phenomenon takes placqfl), we must take the above considerations 
into account. 

From the experiment it is well known that if the radius of the capillary is increased, then the 
superfluidity disappears. Namely, it disappears because collisions with the wall lead to a loss of 
energy. However, the narrower is the capillary, the higher the first transversal mode occurs in a 
narrow capillary. Thus, the energy criterion for which the superfluidity disappears depends on 
the radius of the capillary. 

Therefore, to explain superfluidity in a vessel, we must proceed with the ultrasecond quan- 
tization, and introduce the operators of creation and annihilation of dimers [20j . 



3 Correspondence Principle between the T—p Diagram and the 
Interaction Potential of Lennard— Jones type 

The Van de Waals law of corresponding states establishes the correspondence between differ- 
ent gases. The most accurate calculations can be performed using the methods of molecular 
dynamics, which required a large amount of computer time. Moreover, in these methods, it is 
necessary to know the potential of interaction between particles. 

Here we use the classical scattering problem to establish the law of correspondence between 
the given interaction potential and the diagram p-T, where p is the density and T is the tem- 
perature. This was shown heuristically in [27] and in |32j . 

The isotropy principle (of symmetry in all directions) is one of the key principles of molecular 
physics. It must also be formulated rigorously in mathematical terms as the isotropy principle in 
the theory of Kolmogorov turbulence or the Born-Karman conditions in the theory of crystals 
(the problem of the crystal volume flniteness problem), especially because of rapid development 
of the computer molecular dynamics similar to the computer anisotropic turbulence [7] , [8] . 

The usual argument in molecular physics involves the symmetry of the average motion of the 
molecules in all six directions. Therefore, 1/12 of all particles move toward one another. Since 
there are three axes, it follows that 1/4 of all molecules collide. 

In the two-body scattering problem, the Lennard-Jones interaction potential 

ai2 a6 



^(^''^") = ^' lu. .„||12 - II,. ,„||6 . (28) 



fy*' fy»J 



is usually considered. Here, e is an energy of a depth of the well, a is an effective radius and 
||r' — r"|| is the distance between two particles with radius vectors r' ,r" . In the absence of an 



^In other words, masses of gas move. 



external potential, the two-particle problem reduces to the one-dimensional radial-symmetric 
one. As it is well known [13], in this problem, two quantities (energy E and momentum M) are 
conserved. In the scattering problem, it is convenient to consider, instead of the momentum M, 
another preserved constant, namely, the impact parameter B, so that M = y/EB. 

The scattering of two particles of equal mass may lead sometimes to the formation of a pair. 
If, simultaneously, slight friction or viscosity occurs and a small part of the energy is dissipated 
to radiation (of photons), to noise (of phonons), then, in this formation, it is natural to preserve 
the original preserved (up to infinitely small losses) quantities: the impact parameter B and the 
energy E. 

These heuristic considerations led me to the correspondence principle presented below. 

At each point r E M^, the dressed (or "thermal") potential ip{r) is attractive. In addition, 
because the volume ^ is a large parameter, it follows that if 



(p{r) 




is expanded in terms of 1/V, then 



4 



*(^)=Ci + ^|;£ + 0(F-t), C2>0 (29) 



Since 



11/ // 1 1 2 1 1 / I // 1 1 2 

. , o . , , , o \\v — v\\ \\v -\- v\\ 

\\rf + ||r"||^ = "" / " + "" Y " , (30) 

we can, just as in [2S], separate the variables in the two-particle problem and obtain the scattering 
problem for pairs of particles and the problem of their cooperative motion for r' + r" . 

Then, in the two-body scattering problem, an attractive quadratic potential (inverted 
parabola) is added to the Lennard-Jones interaction potential and the Hamiltonian of the 
scattering problem has the form: 

2 

H = -— + Weff{\\r -r \\), P=\\p-p\\- (31) 

where p' , p" are momenta of two particles with radius vectors r', r" and 

Weffir) = -pr + ^(r), n{r) = As (-^ - -^^ ar\ (32) 

is the effective potential (for short, we denote ||r' — r"|| = r). Here, e is the energy of the well 
depth, a is the effective radius and a = ^. 

Since Hamiltonian B is equal to the total energy, we have 

^--^ = 1^ + ^e//(r) = 1^ + ^^ + u[r\ (33) 

where B is the impact parameter. From ([33]) we derive: 

T) 1X\T] 

^*°*-' = 4m(l - 52/r2) + 1-S7r2- ^^^^ 

Thus, the Hamiltonian splits into two parts: r < B and r > B divided by barrier r = B. One 
of them (for r < B) \s attractive and the other one is repulsive (for r > B). As the repulsive 
barrier, we can consider precisely that noise which creates infinitely small viscosity or friction. 
For r < B, the first term in (I34p is negative while the second term is positive whenever 
a < r < B. The greater the velocity, the less is the energy. The mean velocity is the temperature. 



However, to study the penetration through the barrier of the incident particle, we must plot E 
along the y axis and turn the wells upside down. Then the minimum becomes the barrier and 
the maximum becomes the depth of the well. 

In addition to the attraction problem, there is also the reflection problem for r > B, r < a. 
It is separated by barrier r = B and creates repulsive particles. The repulsive particles obstruct 
the motion of the particles with r < B. This is an analog of the reflection that creates "viscosity" 
in a collisionless plasma. 

In the two-body scattering problem considered above, let us analyse the potential energy 



E{r) = 4e 



,12 



,12 



ar 



B^ 



1 



(35) 



in the attractive case (i.e. for r < B). Here, e is the depth of the well, a is the effective radius 
and B is the impact parameter. By replacing 



B 



B, 



we get rid of a. In what follows, both the wave and the prime will be omitted. 

For a given B and the given potential (j35p . minimum -Bmin ^^(^ maximum -Emax of E{r) are 
achieved at points defined by condition 



^ = 0. 

ar 

For each a < oq, there exists Bo{a) such that -Emax = Emm and, therefore, 

^^ 0. 



dr^ 



(36) 



(37) 



On the graph (a, E) the relation -Emax = Emin determines the analog of the Zeno line. 
Let us now represent the Zeno line on the graph {p,T). 
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Figure 1: T~p diagram for gases corresponding to simple liquids, T^ = T/Tcr, Pr = p/Pcr- The 
Z = j^ = 1.0 line (Zeno-line) on the phase diagram. For states with Z > 1.0 (hard fluids) 
repulsive forces dominate. For states where Z < 1.0 (soft fluids) attractive forces dominate. 
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Figure 2: The curve Z = 1 — ^min for C2 = 1. 



Further, for a fixed p, which is proportional to a, we obtain the asymptotic behavior 



E{r) 



52 



as B 



00, 



where B is the impact parameter, and also the ratio of the difference of the maximal and minimal 
points of E{r) to the maximal point. By our correspondence principle, this ratio corresponds to 
the compressibility factor 

r^ -C/inax -C/min 

~ F 

-'-'max 

and, as -B — 7- 00, for a given p, we obtain the minimum value of Z on the graph {p,T). The 
Zeno line has already been obtained by the rule given above. 

The value of the compressibility factor Z is already plotted along the y axis. We must now 
establish the correspondence with the temperature scale. To do this, consider the ordinate axis, 
i.e., the case p = 0. The point 0.8e, where e is the depth of the well of the Lennard- Jones 
potential corresponds to the Boyle temperature. 

The minimal value of the compressibility factor Z for a given p is equal to 



Zminip) 



E„ 



Er, 



(38) 



B—B„ 



In this problem, eliminating the parameter i?, we found the Zeno line from the condition E-a 
E'mirn and, as S — )■ 00 and for C2 = 1, we see that 



Z = 1 



Er, 



E„ 



(39) 



B— 5>cx) 



This gives a curve shown in Fig. 2 on the plane Z, p or the curve presented in the Fig. 3 in the 
coordinates p/pB- 

We note that the values Z of formula (j39p for i? > 10 remain almost unchanged, i.e., B = 10 
is a large parameter. 

Since for Z < Z^r = 0.29 (the point at which the derivative along the diagonal in Fig. 3 
vanishes), we arrive at a contradiction related to two large parameters B and V, these parameters 
are to be coordinated with one another. 

For C2 = we consider the pair interaction only. The M-tame hyperbola corresponds to this 
"jump". Further near the point Z = 0.444 there must be a short M-wild curve, joined to a new 
M-tame curve, for which C2 = 0. This curve is related to the dimensionless parameter V^^''^B. 

Let us find a point 

Z = 0.444 = ^Z,r. 
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Figure 3: The curve Z = \ — Z^i^ in the coordinates p/pB- 



At this point we have 



if the reduced coordinates 



V = -Vr 



cr: 



2 
T „ P 



T — P 



L J. 



T P 

^cr ^ cr 

are considered. The value of B at the point 1 — Z is equal to 2.271. The critical point Zcr,Pcr 
is connected to this point by a hyperbola, see Fig. 5. 

At the point p = 0.19 there is a modification of the curve of a rather large scale which we 
refer to as the passage from an M — tame curve to an M — wild curve, where M stands for the 
scale. We had expressed a conjecture that this point is related to the so-called triple point. 

After studying the supermodern experimental data, the author resigns the idea that a pres- 
sure jump happens simultaneously to the volume jump. The chaos of clusters which happens in 
the entire conception of the author is not an obstacle to this resignation. The binodals which 
were constructed by the author earlier correspond to the nonequilibrium states and an attempt 
to find a pressure jump with the accuracy up to which thermodynamical problem is solved 
turned out to be futile. 

A thorough analysis of the modern experiments concerning the measurements of the pressure 
of saturated vapor showed that, when at least somewhat receding from the critical point, the 
pressure is equal to the pressure inside the liquid, independently of Earth's attraction which was 
mentioned above. Indeed, at the expense of thermal Brownian motion, when establishing the 
equilibrium at large time, the pressure become equal indeed with sufficiently large accuracy. 

In the table, we present the data corresponding to the resulting diagram (for B = 100 in 
"molecular" units) and note the disagreement between the main dimensionless relations resulting 
from the data from molecular dynamics and the theoretical relations obtained by physicists from 
the the BBGKY hierarchy of equations and the n-particle Gibbs distribution. 

Table 1: 



Zcr 


Pcr/PB 


Tcv/Tb 


0.29 


0.273 


0.39 


0.308 


0.285 


0.38 


0.375 


0.333 


0.296 



The upper row of the table contains the theoretical values of Z^r, Per/ Pb, and T^t/Tb obtained 
by the author, the second row contains the values of these quantities resulting from the newest 
data of molecular dynamics for the Lennard-Jones potential, and the third row gives the values 



based on the Van der Waals equation, which is empirical and even "more" empirical than the 
law of correspondence given in this section. 

For the ratio Tcr/Ts, different data are given in different reference books; this is due to the 
fact that it is very difficult to determine Tb- molecular dynamics does not provide the value 
of Tb- In the physical literature, the interaction potential for inert gases is given; this is the 
Lennard-Jones potential. Therefore, it is interesting to compare the data obtained on the basis 
of the correspondence principle with the corresponding values for gases. 

The value of Zcr was experimentally determined with a high degree of accuracy and it is 
0.29 for inert gases, nitrogen, oxygen, and propane. 

In the correspondence principle given by the author, the value of Pcr/PB (the ratio of the 
critical p to pB, i-e., to the whole length of the interval of p values, where the Zeno line "cuts 
off" the abscissa axis on the graph in Fig. 2 coincides with the corresponding values for water, 
argon, xenon, krypton, ethylene, and a number of other gases. 

Since the physicists (with the exception of V. L. Ginzburg) did not believe the heuristic 
considerations that led to this rule of correspondence, let me present detailed calculations to 
determine the compressibility factor and the analog of the Zeno line for an interaction potential 
U{r) of either of the following types - generalized Lennard-Jones, Buckingham, Kohara, Morse, 
Schommers, Barker. 

The potential energy has the form 

The first derivative is 

^ ('^) - (52 _ ^2)2 • (^^) 

The second derivative is 

E"{r) = —-1—-[2{B^ + W\^)U{r) - 2ar\6B'^ - 3Bh^ + r^) (42) 

+ Ar{B^ - B^r'^)U'{r) + riB"^ - r'^fU''{r)]. 

The first derivative E'{r) is equal to zero if 

_ -2B^Uir) - B'^rU' + r^U\r) 
"- 2r2(r2 - 2^2) ' ^^^ 

The second derivative E"{r) is equal to zero for 

_ 2{B^ + 3B'^r'^)U{r) + 4r{B^ - B^r^)U'{r) + r'^{B'^ - r'^fW (r) 
«- 2r2(654_3^2^2 + ^4), • ( ) 



In order to define the Zeno line, we equate a standing in relation (|43p 

-2B^U{r) - B'^rU' + r^U'{r) 

2r2(r2 - 2^2) 
_ 2{B^ + 3B^r^)U{r) + Ar{B^ - B\^)U'{r) + r2(^2 _ ^2^2jj„ (^^ 
- 2r2(654-352r2 + r4). ' 

This results in the following equation: 

- 8B'^U{r) + 2B'^rU'{r) + r^U'{r) + 2B'^r'^U'' (r) - r'^W (r) = 0. (45) 

Substituting the solution of this equation into relations (jIOj) and ()i3|) . we obtain (in the para- 
metric form) curve E{a) that, in case of the Lennard-Jones potential, is represented by the 
analog of the Zeno line. 



Further, setting E'{r) to be equal to zero and taking B to be sufficiently large {B = 100), 
we get two values of r. Denote by -Emax and -Emim where Smax > ^min, the values of E{r) 
corresponding to these r and consider Z = (.Bmax — -£'min)/-E'max- We determine Tcr{Z, p/pcr) = 
Tcr = 1 in reduced coordinates by using graphs in Fig. 1 for T = Tcr at the point of maximum 
of the binodal Z = Zrr and for Z > Zrr- 
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Figure 4: The binodal and the Zeno line. 



4 Constraint constants 

Van der Waals wrote his remarkable equation of state with two defining constants a and b. 
However, the Lennard-Jones potential describing the interaction between particles (6-12) and 
also containing two defining constants does not lead to the Van der Waals equation either in 
theoretical calculations |27j or in experimental data even for noble (inert) gases. 

Following our ideas based on the analogy with economic laws, we shall construct equations 
of state with reference to three fixed points. 

First, consider the Irving Fisher economic law for assets 



PQ = Mv, 



(46) 



where PQ is the amount of merchandise in its money equivalent, {Q is the merchandise and PQ 
is its money equivalent), M is the money supply, and v is the velocity of circulation. 

Suppose that Mq is some original money supply. (If the Fisher law is considered historically, 
then it is expressed via its gold equivalent.) Let us consider constraints on formula ()46p related to 
the overproduction of merchandise in the form of the following statement. There exist constants a 
and b such that, for 

PQ > bMo or v>c, 

there is no money supply M satisfying relation (I46p . More precisely, the money supply M 
satisfying relation (|46p exists only under two conditions: 



PQ < bMo, 



V < c. 



(47) 



The constraint (I47p is similar to the restriction on the velocity of particles in mechanics. As is 
well known, it led to the revision of the Newton-Galileo mechanics. 
Now consider the equation of an ideal gas 



PV = NT, 



(48) 



where P is the pressure, V is the volume, A^ is the number of particles, and T is the temperature. 
Suppose that Nq is the number of particles as T — t- 0, P — t- 0. Then a similar statement is of 
the form: relation (j48p holds only under the condition 



r>c. 



E = PV> bNo 



where c and b are some constants. 

Relation (j48p for an ideal gas can be represented as 

P = Tp, (49) 

where p = N/V is the density of the gas. Then our constraint can be written as follows: 

For a given gas, there exist constants c and b such that if T > c or P > b, then there is no 

density p satisfying relation (j48|) . 

This question can be stated differently. Suppose that the pressure P and the temperature 

T are given. Does there exist a density p satisfying relation (j49p for a given imperfect gas (or a 

mixture of gases)? 

In modern theory of imperfect gases, it is usual to consider the plane p, T and, in this plane, 

condition ()48p corresponds to the so-called Zeno-line, which, as given by experiment, is a segment 

of the straight line 

p = po(l-^], T<Tb, (50) 



TbJ 

where go, Tb are constants (the constant Tb is called the "Boyle temperature"). 

This dependence was first noticed by Bachinskii, who justified it empirically. Let us present 
the modern general scheme for an imperfect gas in the plane T, p. 

As was already pointed out by the author [31j , the thermodynamic equations of state con- 
taining the quantities: P, the pressure, T, the temperature, p, the chemical potential and, 
respectively, V , the volume, 5, the entropy. A'', the number of particles, constitute a three- 
dimensional Lagrangian manifold in the six-dimensional phase space; moreover, P, T, p play the 
role of coordinates, while V , S, N, respectively, play the role of the corresponding momenta. 

Therefore, for a fixed number of particles A^, the diagram on the plane p = N/V, T given in 
Fig. 1 is the projection on one of the planes of the phase space. Under such a projection, focal 
and caustic points appear. 

The slanting line issuing from the Boyle point in Fig. 1 is called the Zeno line and is of the 
form ([50]) . 

In studying the pressure as a function of the density p, an ambiguity arises on the Zeno line: 

P = .p(l-|). (51) 

Hence the quantity 

p= — + -\/462-46P (52) 

c c 

takes two values for P < b and becomes complex for P > b, a typical simple caustic (or the 
turning point in the one-dimensional quantum Schrodinger equation). 

5 Compression on the plane P, V in Bose— Einstein-type distri- 
butions and the fractal dimensions 

First, let us describe the distribution for Z < Z^r- 

An analog of the potential Vl^ for the number theoretical distribution 7 = (d — l)/2, where 
(i is a "fractal" non-integer dimension and d = D/2, has the form 



Our distribution contains multiplication by a function of V , i.e., the following change occurs 
in the Bose-Einstein distribution: 

V^'^^iy), fj^^X asy^oo. (54) 



The index 7 is sometimes omitted, and hence here it is constant. This function is constant for 
a given dimension. Therefore, Z^i^ obtained for /x = has the form 

where (" is the Riemann function. 

Let us introduce one more quantity. The press P is an intensive quantity conjugate to the 
fractal dimension 7 

P= ' 
For any Z < 1 we have 

^ v^'jv) r(7 + i) r^^^ 



^^{v) r(7 + 2) j;^ 



00 eTde 



e«-«-l 



dj 


















*(k). 


K = 


"r' 


<^7 


(^) 


dip 


(56) 


= 0, we 


obtain 


(F^. 


Further, 


we 


find ;u(y) 


as a 



function of y from the condition Z = 1: 

-^(k) = 1, k = k{V). (57) 



^^^;(^), 



On the other hand, on Z = 1, we have the Bachinskii parabola, which follows from the Zeno 
line. 

P = pbT{1 - ^) (58) 

or 

P = Tbp{1--^), (59) 

PS 

as well as 

T = Tb{1-—). (60) 

PS 

Hence the dependencies P{T), T{p), and P{p) are known: P{T) is the Bachinskii parabola, 
T{p) is a straight line, and P{p) is a parabola. 

Let us find P according to the Bose-Einstein distribution with V replaced by ip^{V), 

p = rr_L — . (61) 

r(7 + 2) Jo e-e^-l ^ ' 

On the Zeno line (i58]l . substituting (|59]l and (f60]) . we obtain the second relation for c^^(F) 

under the additional initial condition that "^ y ' — t- 1 as y — t- 00. The first-order differential 
equation is obtained in parametric form (here k is a parameter). 

The function 'p-yiy) allows us to determine the critical dimension d, and hence 70 = 0.2 as 
follows from ()55p with V = V^r- Therefore, ip^^iy) gives us the solid curve in Fig. 5, which is in 
a good agreement with experimental values for the argon gas component (for P < 1). 

Example 3. The curve Tr = 1 for ideal and imperfect Bose gas in Z, P coordinates. 

The formula for ideal gas is 

—Liyo+i{a) = C(7o + 2), 



where a = a(P) is determined by the equation 

C(7o + 2) 



Li^o+2{a{P)) ^ p 
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Figure 5: The curve Z{p) goes from the point Z = 0.1 to the point Z,,.,. 
curve, after which goes along the thick curve. The hyperbola is dotted. 



0.29 along the thin 



The formula for imperfect gas is 



/(-)Li^„+i(a) = C(7o + 2), 



where a = a{Z, P) is determined by the equation 



P 



The results of computations by the last formula coincide with experimental graph in Fig. 6 for 
Tr = 1 up to the point Z = 0.29. 

How does the curve T^ = 1 go below the point Z = 0.29? 

The compressibility factor Z in this distribution is determined as the following relation 



Z 



PrV 



(62) 



where Pr = P/Pcr, Tr = T/Tcr- Hence, the fractal dimension d is uniquely determined for 



Zj — ^cr — ^cv' 



If we assume that the parameters /3 = 1/T, n, n = /x/T, and 7 are mutually related by the 
condition that the number L of versions of the solution of the Diophantine equations for £ and 
A'^ must increase maximally, then this means that the "specific entropy" S, i.e., the entropy used 
in number theory , takes the maximum values. Since 



7id 



Sid = Z^ — K — )• max. 



yid 

^7+1 






(63) 



where Li is a polylogarithm. 

The condition that there is a relation between /i and 7 completely determines the curve 
Tr = 1 in Fig. 6. The same holds for the curves T^ = const > 1. We write the system of 
equations 

A = dS,d/dfi, 7 = dSid/d-f. (64) 



The limit point is determined by the condition Li^j^i{e'^) = const: 



d_ 
dfj, 



— (Z^+i-k) = 0, 



(65) 




Figure 6: Experimental graph. P^ = Patm/Pcr: ^r = T/T^r, T is the temperature in Kelvin 
degrees, V/R is the volume in crrfi /mole, R is the gas constant, and Z = Pj.V/{RTj.) is the 
compressibility factor. The isochores V/R = const are shown by dotted lines. 



and the constant is determined by the relations on the Zeno line (see Fig. 8). The other boundary 
conditions are determined by the curve in Fig. 3 for Z > Zcr- 

Using the dependence Z{p) {p = \/V) on the thin curve for Z > 0.29 (we mean that, on this 
curve, K is zero) and the dependence on the Zeno line, we obtain that for 7 > 70 the variables 
^'yiy) depend on 7 without a boundary condition of the type ([M|) (for 7 < 70 the variables 
'■P'yiy) do not depend on 7). The boundary condition is adjusted so that the above relations 
and the continuity condition are satisfied. 

Prom (I64p we have 



^7 
dp 



-^ = -TZ, 






7 = 0.2 as p = 0, 



1. 



(66) 



We draw a line between the breakpoint 71 (ju) and the Zeno line at the point Tr = 1 (P 
Pb(1 — I/Tb), Z = l)(jamming effect). The variables (p'^_^ are practically equal to constant. 
By putting '^{p) in 



P 






and 






we obtain the extension of the curve Tr = 1 corresponding to Fig. 6. 

Since it turns out that 7' is close to unity, we present a figure, where 7 linearly depends on V, 
and compare it with the experimental values on the plane Z, P. The character of the behavior 
of 'y{V) and the fluid incompressibility effect as 7 — )■ (the jamming effect) can already be seen 
in this approximation (Fig. 7). 

For Z < 0.17, the jamming effect already takes place, and this permits improve the binodal, 
i.e., the phase equilibrium curve. For 7 < 70, the function 9770 (^) i^ obtained, and hence all 
quantities determined by the potential r?^, can also be determined. 

On the phase equilibrium curve, we know the gas pressure and temperature. They are equal 
to the pressure and temperature of the liquid phase. For Z < 0.17, the picture shown in Fig. 8 
takes place, where p and T do not vary along the rays and are determined on the Zeno line. 
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Figure 7: The isotherms for T^ > 1. The bold hnes denote theoretical isotherms for 7' = 1. 



Since, for the gas branch, we determined the value of T in the triple point, the corresponding 
value of p for the liquid phase we obtain from the equilibrium condition. Thus, we obtain 
ztl^'P'' ^ 0.3 • 10-3, TcV"^'^ ~ 0.55 r,„ p%'P'' ^ 0.7gm/cm3. 

This allows us to introduce one more critical point for the liquid phase and to find (pcr^ ^(^)) 
just as this was done in the case of gas phase, which means that a new potential ^^7^,.^ j^ can be 
introduced in the liquid phase. This is an important fact, useful in liquid dynamics (see [3^). 

Since, in this domain, p and Z < Z^r are related by Z = c/p, where c = const, it follows that 
the decrease in the fractal dimension is determined by the increase in the value of density p, 
and hence also by the increase in the pressure. The constant c is determined by the condition 
c = ZcrPcr- The increase in the pressure can be calculated by the successive approximation 
method. 

In addition, our distribution must be consistent with the Zeno line, which is valid, according 
to experiments, for a wide range of spectra of different gases. Let us show this using an example 
(given below). However, the main thing is that, in our distribution, the dimension is a given 
function of the density p. 

Below the condensate point, for P = 1, the angle of rotation depends on V^/R for V/R = 
Z < Zcr = 0.3. This is a fairly complicated transformation. Nevertheless, since it occurs 
in the plane P, V/R, it does not involve the coordinates T, S/R and preserves the Lagrange 
property, because, in the two-dimensional phase space, any smooth transformation preserves the 
Lagrange property and the coordinates T, S/R remain unchanged. It is important that S is still 
the logarithm of the number of possible variants. 

Let us now pass to the projection of the two-dimensional manifold (obtained above) in the 
six-dimensional space on the four-dimensional phase space. The equality A^ = const cuts a 
two-dimensional surface out of the three-dimensional Lagrangian manifold. The following group 
property holds: as A^ is changed A;-fold, so are the quantities V and S. This implies that this 
manifold is cylindrical. Therefore, it can be projected along p and A^ from the six-dimensional 
phase space on the four-dimensional space, where the explicit action determining the Lagrangian 
manifold is the Gibbs potential $ = Np{P,T). 




Figure 8: Isotherms in the hquid phase region. For Z < 0.17, formula degenerates, and the rays 
T = const stick into the Zeno Hne according to formula ()60p -()63 p . It follows from condition (j63p 
that, for n < 0, the relation between fi and 7 are such that Li^+i(e'^) = const for a fixed T. 
Condition (|65p with (11) for k = N taken into account implies a shift in // of order l/log2 A^- 
This transition on an enlarged scale can be obtained from (I65p and (11). 



Since we have N = const = R, this equality cuts out a two-dimensional cylindrical (by the 
group property) manifold depending on the parameter d = 2j-\-2. Projecting it along N and fi, 
we obtain a two-dimensional manifold in the four-dimensional phase space P, V/R, T, S/R. First, 
consider this projection using Q,lr as an example. 

Here P, T play the role of coordinates, while V/R, S/R that of momenta. Using relation 
P = dQ/dV, we obtain 



P = ^'^T^+'f 



e 



^+^de 



and 



p = if'T^+^ I e^+T 



Q e-A'/^e^-l' 
1 N 



for 7 > 



(67) 



=(C-«^) 



dt 



for 1 < 7 < 0. 



/o i^e^''-"-^ - 1 e(i-'^)^ - 1 

From (j67|) expressing // as a function of P and T: fi = fi{P,T), we obtain the potential 



^ = Rfi{P,T) 



and relations 



V 



9$ 






This gives rise to a potential of the type of the Gibbs potential Rfi{P, T) determining the 
Lagrangian manifold non-uniquely projected onto the plane P, T. 

The usual rule is used for the selection of points at which different branches of the projec- 
tions are joined. The "quantization" of thermodynamics is carried out by the tunnel canonical 
operator, as was described in pO] . 

Just as above, we can project the Lagrangian manifold (after the rotation) in transformed 
(with respect to V and P) coordinates along fi and A^ on the four-dimensional phase space. 
Since the projection does not depend on the order of transformation of the coordinates P and 
V , it follows that the same transformation can be performed after the projection. 

The fact that the potential "mixed" from the free energy and the thermodynamic potential 
arises can easily be explained; see the heuristic example given in |27] . The phase transition 
in liquids is treated there as the creation of "three-dimensional" clusters or domains, as the 
author called them, in which there exists a particle surrounded on all sides by other particles 
(constituting a small coordination sphere). 



Remark 2 The angles of rotation a (in radians) in the plane P, V depending on V are given 
in the following table: 

y > 0.3 a = 0.049 

V = 0.25 a = 0.052 

V = 0.20 a = 0.058 

V = 0.17 a = 0.066 

It should be noted that V = Q.ll corresponds to the last limiting point to which all the values 
V > 0.17 accumulate. This point is a focus; therefore, the "quantization" of thermodynamics 
by the tunnel canonical operator [20j strongly erodes this point. 

Point Z = 0.17 is a focal point, which causes the asymptotics to spread and become smooth. 

Remark 3 In the final formula for the distribution of number-theoretic form, we must also take 
into account the dependence on r that corresponds to the expansion of the thermal potential 
^(r) up to the second power of r in (j29p . 

The phenomenon demonstrated by the graph in Fig. 7 is called jamming. It leads to the 
rapid incompressibility of a fluid, i.e., to the formation of a glass dust, not a crystal (just as in a 
recent eruption of the Icelandic volcano), which is almost insensitive to the pressure, i.e., becomes 
incompressible. For P > 1.5Pcr and at the critical temperature T^ = 1 {Treduced = T/Tcr, and 
hence Z = P/p), line Z = c+dP on the graph (P, Z) implies that p = const, despite the increase 
in pressure. This "new" condensate is not reflected by the diagram T — p m. Fig. 1, where the 
transition to crystal is shown. 

The explanation given to this phenomenon by physicists using the model of hard balls is 
similar to the example of a two-dimensional billiard studied by the author in [26]. Here, as is 
readily verified, the fractal dimension tends to zero along this line, and hence, near this "new 
condensate," it is necessary to use the global asymptotics of number theory presented by the 
author in [26], [32|. 

From our above discussion it follows that for the construction of liquid phase temperature 
curves T^ = const on the graph (P, Z) , one has to use the global asymptotics presented in 
Section 3. We note that, in general, in addition to quantities considered in Section 3, the 
introduced approach allows to determine also the entropy and the fractal dimension. 

The energy spent on the increase in pressure is used not to increase the density, but, rather, 
to effect an internal modification, meaning the decrease in the fractal dimension, and hence the 
decrease in the number of degrees of freedom. In [31], the author called a graph of the type in 
Fig. 7 as a "pit" for the case of economic and revolutionary crises when the revolution in the 
cause of "freedom" leads, as a result of a "new condensate" to the successive decrease in the 
number of degrees of freedom. 

6 A^-particle Gibbs distribution containing all interactions 

Consider a distribution, averaged over different trials, whose number L is much greater than A^, 
of systems of N particles at the same temperature (mean energy), and it is also a distribution 
over energy surfaces H{p,q) = const, p G M^^, q S M^^ (i.e., N = const, L — )• oo). 
This distribution is a distribution over energy levels 

Ei<H{p,q)<E2. 

Let us pass to the statement of the theorem. 

Consider the Weyl quantized self-adjoint operator H = H{p, q) with the discrete positive 
spectrum 

Ai < A2 < • • • < A„ . . . , A„ — )■ 00 as n — ;■ 00 



and the corresponding proper subspaces Pi, P2, . . . , Pn, ■ ■ ■ in the Hilbert space L2(]R^^). We 
assume that the number of eigenvalues less than a given E obeys Weyl's rule as -E — )• 00 (in 
physics, this is expressed as follows: the number of eigenvalues is proportional to the phase 
volume). Sometimes this formula is called Courant's formula. 
Consider the tensor product of L Hilbert spaces 

£2 = ^^2(1^3'^) ® ^2(1^3'^) ® L2iR^^) • • • ® L2(M=^^), (69) 



where L — t- 00, and define the operator 

Hl = # ig) 1 (g) 1 • • • 1^ 1 + 1 (g) ^ 1^ 1 • • • (8) 1 + • • • + 1 1^ 1 • • • 1 (8) ^^ . (70) 

L 

Determine the parameter Be from condition 

Sp( He-^E^) 



Sp(e 



-bEH\ 



E. (71) 



Denote Lq = Sp(e ^ ). The occupation numbers of the operator H^ corresponding to 
the eigenfunction ^j of the operator H less than or equal to EL, are denoted by Lj. The 
sum Li + Lj_|_i + • • • + -Lj+fc is denoted by Ll , the occupation numberqj corresponding to the 
collection of functions ^j, ^j+i, . . . , ^j+n, while the projection operator on the subspace spanned 

(k) 

by proper subspaces corresponding to them is denoted by P^ . 
For a Gibbs ensemble, the following theorem holds [24] . 

Theorem 1 There exist constants Ci, I = 1,2,..., such that, for all i > 1 and n > 0, the 
occupation numbers L^"' satisfy the inequality 



4"^ - BSpiPy^'^e-'^'') > B^/LolnLoi^iLo)) < QL^', 1 = 1,2,..., (72) 



where B = L/Lq, ip{x) is a positive function tending arbitrarily slowly to +00 as x ^ +00, and 
P(-) is the ratio of the number of eigenvalues of the Hamiltonian Hl which do not exceed LE and 
satisfy the inequality in the parentheses to the total number of eigenvalues of the Hamiltonian 
Hl which do not exceed LE (the spectrum density]^. 

Since ^ XiLi is the eigenvalue of the operator Hl corresponding to the eigenfunctions, it 
follows that the sum ^ Lj over all Li is the number of all eigenvalues of the operator Hl less 
than LE. 

In) 

Theorem 1 implies that if we prescribe L^ on the interval Aj , . . . , \i+n and calculate the 
ratio of the number of eigenvalues with given Li, Lj+i, . . . , Lj_|_„ and other arbitrary Lj for 
j<i, j>i + n to the total number of eigenvalues of the Hamiltonian Hl not exceeding LE, 
then this ratio tends to zero as Lq '^ for any s outside the interval 



Lf ^ ~ SpiPre-'^"") ± B/LolnLo^(Lo) 

for Theorem 1. 

Here we are referring to the spectrum density of the Gibbs ensemble and to the fact that for- 
mula (j7ip determines the mean value of the energy and the value of LE bounding the spectrum 
density of the operator Hl. We can also say that this quantity is the total number of eigen- 
values of the operator HlPle-, where Ple is the projection operator onto all proper subspaces 



To define the number L^ \ whose physical meaning is the number of particles, we can assume that the 
particles are indistinguishable. 

^Landau and Lifshits used the term "spectrum denseness" [T^], P- 44 for the microcanonical distribution. 



corresponding to all Hn < LE, where the /in are all the eigenvalues of the operator H^ not 
exceeding LE. 

Suppose that V^ is the projection operator onto all proper subspaces of the operator PleHl 
lying outside the interval enclosing the eigenvalues of the operator 



Then the ratio of the number of eigenvalues of the operator V^Hl to the number of eigenvalues 
of the operator PenHl tends to zero faster than any power of Lq . 

To pass to the classical Gibbs distribution, it is not necessary to pass to the limit as /i — ?• 0. 
It suffices to define the self-adjoint operator with eigenvalues equal to those intervals of the 
phase volume that are satisfied by the eigenvalues. Then we only have to pass from the sums 
to integrals by using the Euler-Maclaurin formulas. 

Let us split the phase space {p, q) E M.^^ into a finite number of domains 

Ei<H{p,q)<Ei+i, (73) 

where I = 0, . . . ,s - 1, Eq = 0, Es = E, p G R^^ , q £ M^^, and M^A^i ig the phase space with 
coordinates Pi,qi, P21Q2, ■ ■ ■ , PNi Qn- Let us perform an ordered sampling with replacement of 
Li from the partition of the domains of the space M^^^ into the "box" Ei < H{p,q) < -E^+i, 
under the condition 

L 
Y,H{puqi)<LE. (74) 

j=i 
Prom the physical point of view, ordered sampling means that we consider L distinguishable 
SAT-dimensional particles. Let p^ denote the cluster "density" in the energy interval Ei < 
H{p,q)<Ei+^. 

Suppose that the above conditions on the function H{p, q) hold. Let L be given. Let us 
determine b from the condition 

00 -| 

^-bH{p,q) ^p^q = i^ b=--. (75) 

^-^ 

We define E in ([71]) as 

POO 

/ H{p,q)e-^"^P^'^Updq. 
Jo 

Then the following theorem is valid. 
Theorem 2 The following relation holds: 



• Ei<H{p,q)<Ei^-l 

where m is any integer. 



I e-^^^P''''>dpdq >VLhn^ij{L)] KL-"", (76) 



Here the probability P is the Lebesgue measure of the phase volume given in parentheses 
76l) with respect to the whole phase volume (f74j) . 

The proof of Theorem 2 is given in [29], Section 2; see also 



7 On the Maxwell— Boltzmann distribution 

The distribution of the Bose-Einstein type distribution is determined from the relation 

£, peM^ qe 



H{p,q)dpdq _ 3 3 



& 



l3{Hip,q)+K) _ I 



dpdq = —, (77) 



^f(H{p,q)+K) _l\ ^ ^ y 



where 



1 u v 



If ;U — )• — oo, we get the Maxwell-Boltzman distribution e~^P /2m+u{q))/T _ 

On the assumption that the potential field varies very slowly, the function of the coordinates 
has the form 



'V, 

Only on this assumption, we can pass from the Maxwell-Boltzmann distribution e~'^ /2m+u{q))/T 
(by integrating it over the momenta) to the Boltzmann distribution of the form 

e-/3«(<?)^ (78) 

because the Maxwell-Boltzmann distribution is not a distribution of the particle number density 
with respect to momenta and coordinates. It only gives the number of particles between the 
energy levels 

h u{q) = const. 

2m 



Therefore, only if u{q) is of the form 



<Q) = ^ > 3/^ 



then this can be done in the thermodynamic limit. 

Here we must not forget that the Boltzmann distribution is not also a distribution with 
respect to coordinates q, but rather with respect to the level surfaces of function u{q). 



8 Conclusions 

We began with usual statistics used in molecular physics (see Section 3 before formula (24)). 
We determined Zcr and Zeno line. But only Zeno line {Z = 1) and the critical points are stable 
under transition to statistics of identical particles (monomers). The critical points give a very 
strong focus, and Z = 1 means the "complete" victory of Shannon's entropy. We obtain a noble 
gas out of glassy dust (the hydrodynamics of such a dust was studied in ^20j ) when most of its 
particles are distinguishable from each other. Namely, in noble gases, there are no preferences 
for dimers, trimers, etc. Therefore, thermodynamics of noble gases gives good models of the 
mixture of Bose-Einstein statistics and Boltzmann statistics. There is a remarkably simple law 
of the mixture of statistics here. As Z = 1, this law is called the Zeno line. 

In the construction in which we define the critical temperature, we consider the difference 
of the energies of the stable and unstable rest points. Furthermore, we use only small friction 
and viscosity. Indeed, as a result of small viscosity, a classical particle, after having flown "just" 
above the barrier, will lose its energy and, on reflection, will hit the barrier and will continue 
hitting the walls of the barrier and the walls of the potential in the form of a well until it 
precipitates to the bottom. (The value of the energy of the stable rest point is Ei and the value 
of the energy of the unstable rest point E2.) To knock out this particle from the well trap, 
the required kinetic energy must be "slightly" greater than the difference between the energies 
El — E2 of the rest points. This can interpreted as a break-up of the dimer in the collision with 
a fast monomer. The equilibrium is violated when the depth Ei — E2 of the well decreases and 
its width increases as the absolute value of the energy Ei of the stable rest point decreases (as 
the impact parameter increases). It is natural to regard the increasing width as a result of an 
increase in the numbers of dimers in the trap (in quantum theory, this corresponds to an increase 
in spectral density). Dimers clusters can survive and equilibrium between the monomers and 



the dimers can be preserved only if the dimers themselves combine into cluster domains and 
create their own barrier, which is a microanalog of a surface film. 

In problem (30) as C2 = 0, different barriers and wells occur for different values of B. At the 
rest points, -Emm and -Bmax! the velocity is zero; therefore, they can be determined only from 
the potential term. 

We are now dealing not with just one particle, but with a pair of particles whose mass center 
is in the trap. Therefore, the difference -Bmax — -E-min is the energy needed to knock out this pair 
(the dimer) from the trap. 

Experimentally, we can calculate the percentage of dimers in a gas. It is clear how dimers 
are created and split by monomers. Further, their mean number is calculated. The higher is 
the temperature, the greater is the mean energy of monomers and the smaller the number of 
dimers. 

By our calculations, we have Tct/Tb = 2.79. According to the contemporary handbooks, 
this value is Tcr/Ts = 2.72 for argon (Ar), Tcr/Ts = 2.71 for krypton (Kr), Tcr/Ts = 3.157 for 
methane (CH/^), and Tcr/Ts = 2.6 for nitrogen (A^2)- (Other data is given in pj.) 

In Table below, comparative data for T^r/^ are given. 



Substance 


e,K 


T,./4 


Ecr ■ e/k 


Ne 


36.3 


11 


10.5 


Ar 


119.3 


37 


35 


Kr 


171 


52 


50 


N2 


95, 9 


31 


28 


CHi 


148.2 


47 


43 


C2H6 


243.0 


76 


70 



Since, as is well known, Z begins on the {p, T) coordinates, we can use the dependence (given 
above) up to Z = 0.444 and then include the thermic potential, because we must also take into 
account the influence of a third particle. As was already stated in previous papers, the dressed, 
or "thermic," potential ^(r) is attractive. 

The "mixture" of statistics allowed us to solve the famous problem of Gibbs paradox as 
a contradictory example in the system of phenomenological axioms of thermodynamics. As 
Poincare noted, we turned to "arithmetics" to solve this problem and used the mixture of the 
statistic from the number theory and the usual statistic. In numerous attempts to solve this 
problem were base on the use of a mixture of particles [5] and their internal structure. The fact 
that there are foci at the critical point and the Zeno line is stable (i.e., the compressibility factor 
is equal to unity) allowed us to do this. 

This problem was attacked not only by many physicists but also by philosophers, for the 
first time, by Poincare in his philosophic works, and also by B. V. Kedrov, S. D. Haitun, and 
many other. 

I wish to thank D. S. Minenkov and A. V. Churkin who performed the computations. I also 
wish to express deep gratitude to Professor V. S. Vorobiev for his enthusiasm, for very productive 
discussions, and for his efforts to verify the computational results and to compare them with 
experimental data. I am cordially grateful to philosopher academician A. A. Guseinov for very 
useful consultations. 
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